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Abstract. The enhanced binary tree (EBT) is a nontransitive graph which has two 
percolation thresholds p c \ and p C 2 with p c \ < p C 2- Our Monte Carlo study implies 
that the second threshold p C 2 is significantly lower than a recent claim by Nogawa and 
Hasegawa (J. Phys. A: Math. Theor. 42 (2009) 145001). This means that p c2 for the 
EBT does not obey the duality relation for the thresholds of dual graphs p C 2 + P C \ = 1 
which is a property of a transitive, nonamenable, planar graph with one end. As in 
regular hyperbolic lattices, this relation instead becomes an inequality p C 2 +p cl < 1. 
We also find that the critical behavior is well described by the scaling form previously 
found for regular hyperbolic lattices. 
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Figure 1. (a) Mass of the root cluster, so, divided by L, the number of generations 
in the EBT. The crossing point indicates the first percolation transition point, p c i, 
where an unbounded cluster emerges. Inset: Scaling collapse by equation ([1]) with 
Pd = 0.304 found in [1]. (b) The number of boundary points connected to the root 
node, denoted as b, also shows a crossing point at p = p c \. Inset: Scaling collapse by 
equation ([2]) with the same p c \ as above. 



Recently, Nogawa and Hasegawa pQ reported the two-stage percolation transition 
on a nonamenable graph which they called the enhanced binary tree (EBT). While the 
first transition had little ambiguity, they mentioned that the behavior at the second 
threshold did not look like a usual continuous phase transition. 

A quantity of interest was the mass of the root cluster, denoted as so, where the 
root cluster was defined as the one including the root node of the EBT. Using this 
observable, we briefly check the first transition point, p&, where an unbounded cluster 
begins to form. As in pj, we have used the Newman-Ziff algorithm [3l H] and taken 
averages over 10 6 samples throughout this work. The number of generations, L, of the 
EBT defines a typical length scale of the system, and p] showed the finite-size scaling 
of s as 

8 /L<xf 1 [(p-P*)L 1/v ], (1) 

with v = 1. Figure [T]( a) confirms both of the percolation threshold p c \ and the scaling 
form, equation (1T1). Equivalently, one can measure b, the number of boundary points 
connected to the root node, which becomes finite above p c \ as shown in figure [H( b). It 
also scales as 

b<xf 2 [(p-Pd)L 1/v ], (2) 

with the same exponent v. Comparing this with [2], we see that the percolation 
transition in the EBT at p — p c \ belongs to the same universality class as that of regular 
hyperbolic lattices. One may argue that this scaling form actually corresponds to the 
case of Cayley trees [2]. The convincing results in figure [U imply that the estimation in 
|T] for the dual of the EBT, p cl = 0.436, is also correct. 

On the other hand, the second percolation transition at p — p C 2 indicates uniqueness 
of the unbounded cluster. We have thus employed a direct observable to detect this 
transition, i.e., the ratio between the first and second largest cluster masses [2]. The 
idea is that even the second largest cluster would become negligible if there can exist 
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Figure 2. (a) Ratio between the second largest cluster mass, S2, and the first largest 
one, s\. The crossing point lies at p w 0.48, which implies that only one cluster will 
remain dominant in the system. ( b) Number of boundary points connected to the root 
node, b, divided by the total number of boundary points, B. The dotted black curve 
marked by 00 indicates the extrapolation result from equation ([3]). (c) Scaling collapse 
according to equation (jlj), where we set p C 2 = 0.48, <j) = 0.84, and \jv = 0.12. 



only one unique unbounded cluster. Measuring s%/s\ in the EBT, where Sj means the 
zth largest cluster mass, we have found the second transition at p C 2 ~ 0.48 (figure a)), 
certainly lower than Nogawa and Hasegawa's estimation, p = 0.564. 

As an alternative quantity for p&, we divide b by the number of all the boundary 
points, B. This fraction b/B is supposed to become finite above p C 2 [2]. Based on the 
Cayley tree result [2], we have assumed that as the system size N varies, one can write 
down the following asymptotic form: 

b/B ~ Cl N^ + pa, (3) 

with some constants c\ and C2 and an exponent <p. From the finite-size data, 
we extrapolate the large-system limit by equation (ED, which suggests p C 2 ~ 0.49 
(figure Wyb)). This is very close to the estimation above from S2/ S\. Moreover, in 
accordance with equation (jHJ), we have suggested the following scaling hypothesis to 
describe the critical behavior at this transition point [2]: 

b/B<xN*- 1 U{p-Pc2)N l /% (4) 

with an exponent v. Applying this hypothesis to EBT data, we see that <fi = 0.84 and 
1/p — 0.12 give a good fit (figure E](c)) with the same value of p C 2 = 0.48, where the 
numeric values of the scaling exponents are again consistent with [2]. 

To make a direct comparison to the observation in PQ, we have also calculated 
the mass fraction of the root cluster, Sq/N, as a function of p. As above, performing 
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Figure 3. Mass fraction of the root cluster. The extrapolation result is represented 
by the dotted black curve named as oo. The arrow indicates p = 0.564, predicted as 
the transition point in pQ. 

extrapolation to the large-system limit, we see that this quantity becomes positive finite 
at p > 0.49 (figure EJ. 

All of these observations suggest that the predicted value of p c2 in [1] is too high, 
and it seems that this overestimation led them to consider 'discontinuity' since sq/N 
became already so large at that point as shown in figure O 

Finally, even though our estimation suggests such a different p c2 that p C 2 + p c i < 1, 
we note that it does not violate the duality relation proved in [3] for a transitive, 
nonamenable, planar graph with one end: As Nogawa and Hasegawa correctly pointed 
out P, the EBT does not possess transitivity. The inequality p C 2+P c i < 1 was explicitly 
verified for a pair of hyperbolic dual lattices {7,3} and {3,7} in [2]. This inequality 
means the existence of a narrow region of p between p c2 and 1 — p cl , where one would 
find a unique unbounded cluster in a given graph whereas infinitely many unbounded 
clusters in its dual graph. Such a region does not exist for a transitive case [5]. A 
typical state in this region is illustrated in figure HI which shows a situation with many 
unbounded clusters of radii comparable to L at the same time as a single unbounded 
cluster occupies the dominant part of the dual graph. 
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Figure 4. Visualization of a triangular hyperbolic lattice projected on the Poincare 
disk, where the maximum length from the origin is chosen to be L = 4. Bonds are 
randomly occupied with probability of p — 0.42, which are colored red, while only the 
rest of them appear as occupied in the dual lattice, as colored green, so that the dual 
probability corresponds to p = 1 — p = 0.58. Note that p lies between p C 2 and 1 — p cl , 
since this structure has p C 2 ~ 0.37 and its dual has p cl « 0.53, according to [2]. While 
most clusters have been already absorbed into the largest red cluster, many of green 
clusters still have radii comparable to L since p cl <p< p c2 ~ 0.72. 



